#pragma once

#include "para.h"
#include "force.h"
#include "mol_sys.h"

void Compute_long_range_corrections(double& long_range_correction_energy, double& long_range_correction_pressure);
void inside_box(std::vector<Atom>& mol);
void remove_com_velocity(std::vector<Atom>& mol);
void rescale_velocities(std::vector<Atom>& mol);

/**************************************************************************PIMD**************************************************************************/
void inside_box_PIMD(std::vector< std::vector<Atom> >& mirr);
void remove_com_velocity_PIMD(std::vector< std::vector<Atom> >& mirr);
void rescale_velocities_PIMD(std::vector< std::vector<Atom> >& mirr);